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Abstract This paper considers volume averaging in the quasispherical Szek- 
eres model. The volume averaging became of considerable interest after it was 
shown that the volume acceleration calculated within the averaging frame- 
work can be positive even when the local expansion rate decelerates. This 
issue was intensively studied within spherically symmetric models. However, 
since our Universe is not spherically symmetric similar analysis is needed in 
non-symmetrical models. This papers presents the averaging analysis within 
the quasispherical Szekeres model which is a non-symmetrical generalisation 
of the spherically symmetric Lemaitre-Tolman family of models. In the qua- 
sispherical Szekeres model the distribution of mass over a surface of constant 
t and r has the form of a mass-dipole superposed on a monopole. This paper 
shows that when calculating the volume acceleration, d, within the Szekeres 
model, the dipole does not contribute to the final result, hence d only de- 
pends on a monopole configuration. Thus, the volume averaging within the 
Szekeres model leads to literally the same solutions as those obtained within 
the Lemaitre-Tolman model. 

PACS 98.80-k, 95.36.+X, 98.65.Dx 



1 Introduction 

In the standard approach to cosmology it is assumed that the Universe can be 
described by the homogeneous Friedmann model. Within such a framework in 
order to correctly describe cosmological observations one needs to postulate 
the existence of dark energy, which in its simplest form can be represented 
by a cosmological constant. However, our Universe on scales up to at least 
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100 Mpc is very inhomogeneous. Thus, it can evolve differently than the 
homogeneous model. The difference between evolution of the homogeneous 
model and the inhomogeneous Universe is known as a backreaction effect. 
The direct study of the dynamical effects of inhomogeneities is difficult since 
both general matter distribution and the numerical evolution of cosmological 
models employing the full Einstein equations are unavailable at the level of 
detail which would make them useful in studying this problem. Currently, 
therefore, one of the most popular approaches to backreaction is via averaging 
methods. 

In the averaging approach to backreaction, one considers a solution to the 
Einstein equations for a general matter distribution and then an average of 
various observable quantities is taken. If several assumptions are introduced 
(see Sec. the averaging leads to the Buchert equations. For a review on 
backreaction and the Buchert averaging scheme the reader is referred to O 
U|. Within the Buchert averaging scheme there are a number of examples 
where it was shown by using spherically symmetric inhomogeneous models 
that one can obtain negative values of the volume deceleration parameter 

However, we should be aware that the results of averaging can be gauge- 
depended Another problem regarding an application of spherical sym- 
metric models is the problem with the age of the Universe. Within the models 
studied in [XOj those of realistic density distribution and with q < had very 
large values of the ts function, of amplitude 10 10 y (this means that the 
age of the Universe within such models is unrealistically smalQ). This fea- 
ture, however, can be an artefact of assumed spherical symmetry. Therefore, 
it is of great importance to study averaging in non-symmetrical inhomoge- 
neous models. One of the immediate candidates is the Szekeres model [17] . 
The Szekeres model is a generalisation of the Lemaitre-Tolman model that 
has no symmetry [20] . Within the quasispherical Szekeres model one can de- 
scribe two [21] or even three structures [22] , Thus, the Szekeres model not 
only allows us to study how cosmic structures affect their evolution but also 
enables the examination of the volume acceleration of such systems which 
consist of several structures. This paper, therefore, addresses the subject of 
volume averaging in the Szekeres model. 



2 Buchert equations 

Since the Buchert averaging scheme involves the volume average it applies to 
averaging of scalars only - volume averaging of tensors leads to noncovariant 
quantities (for a review and a detailed discussion about tensor averaging the 
reader is referred to [IS]). The average of a scalar & is then equal to 




(1) 



1 This, however, does not apply to the two-scale averaging ap proach. Fo r details 
on the two-scale models see works by Rasanen [12] or Wiltshire [T3|I14II15| . 
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where h is a determinant of a 3D spatial metric, h a p = g a p — u a up; I? is a 
domain of averaging, and Vd its volume, 

V v = [ d 3 xy/^h. (2) 
Jv 

The Buchert equations are obtained if the following are assumed 

— the Universe is filled with an irrotational dust only, 

— the metric is of the form ds 2 = dt 2 — gijdx l dx^ (3+1 ADM space-time 
foliation with a constant lapse and a vanishing shift vector). 

Then by averaging the Raychaudhuri equation we obtain [I] 

3— = -Mp)v + Qv, (3) 

where a dot (') denotes dt ; the scale factor av is defined as a cube root of the 
volume 

av = (Vv/V Vz ) 1/3 , (4) 
(where Vx> ; is an initial volume) ; and the backreaction term Qv is given by 

Q v ^l({0 2 ) v -{0) 2 v )-2{a 2 )v, (5) 

where is the scalar of expansion and a is the shear scalar. Averaging of the 
Hamiltonian constraint leads to [I] 

3?f=8n{p) v -hK) v -±Q v , (6) 

where (TZ)v is an average of the spatial Ricci scalar ^IZ. The above is 
compatible with ((3j) if the integrability condition holds 

^fdt (Qva v ) + 4-3* ((R)va 2 v ) = 0. (7) 

a T> a T> 

Equations ([3]) and ((6j) are very similar to the Friedmann equations, where 
Qt> = 0, and p and 1Z depend on time only. In fact, it can be shown that 
they are kinematically equivalent to Friedmann equations with a scalar field 

Using ([3]) and iJU we can calculate the deceleration parameter 

= dyay -47rG(p) p + Q-p 

q ~ d% ~ 87rG{p)v-k(K)v-hQv' U 
In the standard approach to cosmology where Friedmann models are em- 
ployed the case of q < implies that A > 0. However, as it was shown within 
inhomogeneous but isotropic models [5,6,7,8,9,10] we can have q < even 
if A = 0. This suggest that the apparent acceleration of the Universe might 
be explained not by invoking dark energy but by taking into account mat- 
ter inhomogeneities. In Sec. [4] we will examine this issue by employing the 
non-symmetrical Szekeres model. 
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3 The Szekeres model 

The metric of the Szekeres model is of the following form |17| 

ds 2 = dt 2 -X 2 dr 2 - Y 2 {dx 2 + dy 2 ). (9) 

For our purpose it is more convenient to adopt a pair of complex conjugate 
coordinates 

C = x + iy, ( = x-iy, (10) 
so that the metric becomes 



where 



and 



ds 2 = dt 2 - X z dr 2 - F^dCdC. (11) 



£ = a(r)(( + b{r)C + c(r)( + d(r), e = 0, ±1. 



Here a prime (') denotes d r . 

The case where e = — 1 is often called the quasihyperbolic Szekeres model 
(for a detailed discussion on the quasihyperbolic Szekeres models see [23J), 
e = quasiplane (for details see [231124] ) . and e = 1 quasispherical (for details 
see [25]). Although it is possible to have within one model quasispherical and 
quasihyperbolic regions separated by the quasiplane region [23J, only the qua- 
sispherical case will be considered here. This is because the averaging within 
the quasihyperbolic and quasiplane requires a special treatment. Firstly, an 
area of a surface of constant t and r in the quasihyperbolic and quasiplane 
models is infinite. Secondly, there is no origin - in the quasihyperbolic model 
r cannot be equal to 0, and in the quasiplane r can only asymptotically 
approach the origin, r — * [23] . 

In the quasispherical Szekeres model a surface of constant t and r is a 
sphere of radius #(r, t) [17]. Thus, the quasispherical Szekeres model is a 
generalisation of the Lemaitre-Tolman model [18, 19J. Within the Szekeres 
model shells of matter, however, are not concentric. The quasispherical Szek- 
eres model becomes the Lemaitre-Tolman (LT) model when £' = 0. In this 
case 

X = —=—, Y = -, JlA^d0 2 + siii 2 

Originally, Szekeres considered only a case of p = = A. This result 
was generalised by Szafron [26] to the case of uniform pressure, p = p(t). A 
special case of this solution, the cosmological constant, was discussed in detail 
by Barrow and Stein-Schabes [27]. In the case of p = = A, the Einstein 
equations reduce to 

*> = 2 -f-*, (12) 
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M' - UI£'/£ 
A *P = &(*>-*£>/£) > (13) 

where M(r) is another arbitrary function. In a Newtonian limit M is equal 
to the mass inside the shell of radial coordinate r. Although the p func- 
tion, as seen from (fT3|) . is a function of all coordinates, it can be shown 
that the distribution of mass over each single sphere of constant t and r 
has the structure of a mass dipole superposed on a monopole, p(t, r, f , Q = 
Pmon(t,r) + pdi P (t 7 r 7 (,C) |28,29,30j. In general case, the orientation of the 
dipole axis is different on every constant- (t, r) sphere. 

The dipole contribution vanishes when £' = and then the Szekeres 
model reduces to the LT model. 

The scalar of expansion is equal to 

_ ^+2gg7g-3gg7g 
U ~ U ;Q _ &- $£>/£ • [M) 

The scalar of shear is 

2 



a 2 = - 



3 y& -$£'/£ J ' 
The spatial Ricci scalar ^TZ is equal to 

(») V -n k ( 4>k'/k-2$£'/£ 

In the LT limit these scalars reduce to 



(15) 



(16) 



^ & & , i / & S\ ,„ ($ky , s 



4 Averaging in the quasispherical Szekeres model 

This section considers the volume averaging within the quasispherical Szek- 
eres model. The volume is calculated around the observer located at the ori- 
gin. It will be shown that the dipole configuration does not affect such quan- 
tities as (p)x>, Qt>, or (lZ)v- These functions only depend on the monopole 
configuration. First let us notice that the volume in the Szekeres model is 
exactly the same as in the LT model (i.e. if £' = 0). Following the method 
presented in [28] we obtain 



Vd 



dr 



dCdC XY 



dr 



dcdc _£l_( # ._4)_L 



dr 



dCdC 1 <P 3 d 

~£^ + 2^T = %dr~ 



dCdC\ 

£ 2 J 



(18) 



c 



Since d(d(/£ 2 is the metric of a unit sphere, hence 

dCdC 



£ 2 



Air. 



Thus 



V v = Att I dr (19) 

J Vl — k 

o 

The same result is obtained if initially £' is set to zero. Thus, the dipole 
component does not contribute to total volume. As it will be shown below it 
also does not contribute to (p)x>. 

r-D 







1 f dr f f *cd(^_ (V _ $£\ M ' - 3M£l i £ 



4ttVd J J J £ 2 Vl~^k V £ J <P 2 (<P' - $£'/£) 

o 

Vd J VT^k 8ttV v J VT^k dr \J J £ 2 J 

o o 

1 f M 1 

o 

The average of the scalar of the expansion is 

j_ 7 r r dCdC ( £\ *> + 2fi*>/*-3*e>/e 

[ >v VdJ J J £ 2 VT=k\ £J &-*e'/e 





ru „ / . . \ 

= / dr^= — + 2- . (21) 

o v 7 

As above, £' /£ does not contribute to the final result and, as seen from $T7\ . 
the result is the same as in the LT model. To calculate the backreaction term 
Qt> we still need to find the average of the following quantity 

dCdC 'I' 1 (',,, 







- <p— 

£ 2 VT^k, V £ 



<P' + 2<P<P'/<P - W£' '/£) - (&' - <P<P'/<P 



{& - $£'/£f 
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2 1 f , f /"dCdC <P 2 L$ 2 - n $ 2 ^£' n $ 2 £'\ 

o x 7 

v 7 

This result is also the same if initially £' is set to zero. Finally let us notice 
that (1Z)t> is also the same as in the case of £' — 0. Averaging relation lfl6|) 
yields 







I'D 



' dr ) 1 ■ (23) 







Inserting l|20[l - lj23"|) into (J3j) and (J6j) we see that the dipole configuration 
contributes neither to a-p nor to d-p. The volume acceleration, as well as the 
volume deceleration parameter, depends only on the monopole distribution 
— it depends only on functions of t and r [i.e. M(r), k(r) and <P(t,r)], and 
not on £ and C [the dependence on £(t,r, £, £)) vanishes]. Moreover, (p)x>, 
Qx>, and (7?.)x) are literally of the same form as if initially £' is set to zero 
(the LT case). 



5 Conclusions 

In this paper the volume averaging within the quasispherical Szekeres model 
has been investigated. The Szekeres model is a generalisation of the LT model. 
The density distribution in the quasispherical Szekeres model has the struc- 
ture of a time-dependent mass dipole superposed on a monopole. When cal- 
culating the volume acceleration (a to be more exact) or volume deceleration 
parameter (q) within the quasispherical Szekeres, the dipole does not con- 
tribute to the final result and a only depends on a monopole configuration. 
The solutions are the same if initially £' was set to zero, thus the results 
and conclusions found when studying averaging within the LT models also 
apply to the Szekeres models. For example, we can, without any further cal- 
culations, conclude that Qx> = when fc = 0. This is an implication of the 
result obtained by Paranjape and Singh [6] who showed that in the parabolic 
(fc = 0) LT model the backreaction term, Qx>, vanishes. Another result ob- 
tained within the LT models which, as has been shown, also appears to apply 
to the Szekeres models is that Qx> > is only possible for unbounded sys- 
tems, k < [2]. However, as shown in [10], in most cases this also requires 
that the bang time function, is, is of amplitude of 10 10 y. 
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